Beyond Tonks-Girardeau: strongly correlated regime in quasi-one-dimensional Bose 

gases 
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We consider a homogeneous ID Bose gas with contact interactions and large attractive coupling 
constant. This system can be realized in tight waveguides by exploiting a confinement induced 
resonance of the effective ID scattering amplitude. By using a variational ansatz for the many-body 
wavefunction, we show that for small densities the gas-like state is stable and the corresponding 
equation of state is well described by a gas of hard rods. By calculating the compressibility of the 
system, we provide an estimate of the critical density at which the gas-like state becomes unstable 
against cluster formation. Within the hard-rod model we calculate the one-body density matrix 
and the static structure factor of the gas. The results show that in this regime the system is more 
strongly correlated than a Tonks-Girardeau gas. The frequency of the lowest breathing mode for 
harmonically trapped systems is also discussed as a function of the interaction strength. 
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The study of quasi- ID Bose gases in the quantum- 
degenerate regime has become a very active area of re- 
search. The role of correlations and of quantum fluctua- 
tions is greatly enhanced by the reduced dimensionality 
and ID quantum gases constitute well suited systems to 
study beyond mean-field effects Q]. Among these, par- 
ticularly intriguing is the fermionization of a ID Bose gas 
in the strongly repulsive Tonks-Girardeau (TG) regime, 
where the system behaves as if it consisted of noninteract- 
ing spinless fermions 2] . The Bose-Fermi mapping of the 
TG gas is a peculiar aspect of the universal low-energy 
properties which are exhibited by bosonic and fermionic 
gapless ID quantum systems and are described by the 
Luttinger liquid model |3J • The concept of Luttinger liq- 
uid plays a central role in condensed matter physics and 
the prospect of a clean testing for its physical implica- 
tions using ultracold gases confined in highly elongated 
traps is fascinating 
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Bosonic gases in ID configurations have been real- 
ized experimentally. Complete freezing of the transverse 
degrees of freedom and fully ID kinematics has been 
reached for systems prepared in a deep 2D optical lat- 
tice 0, Q. The strongly interacting regime has been 
achieved by adding a longitudinal periodic potential and 
the transition from a ID superfluid to a Mott insulator 
has been observed A different technique to increase 
the strength of the interactions, which is largely employed 
in both bosonic and fermionic 3D systems [!j but has 
not yet been applied to ID configurations, consists in 
the use of a Feshbach resonance. With this method one 
can tune the effective ID coupling constant gm to es- 
sentially any value, includin g zfc oo, by exploiting a con- 
finement induced resonance [lfj. For large and positive 
values of giu the system is a TG gas of point-like im- 
penetrable bosons. On the contrary, if g\£> is large and 
negative, we will show that a new gas-like regime is en- 



tered (super- Tonks) where the hard-core repulsion be- 
tween particles becomes of finite range and correlations 
are stronger than in the TG regime. Some consequences 
of this new regime on the energetics of small systems in 
harmonic traps have already been pointed out in a pre- 
ceding study In this Letter we investigate using 
Variational Monte Carlo (VMC) techniques the equa- 
tion of state and the correlation functions of a homo- 
geneous ID Bose gas in the super- Tonks regime. We find 
that the particle-particle correlations decay faster than 
in the TG gas and that the static structure factor ex- 
hibits a pronounced peak. The momentum distribution 
and the structure factor of the gas are directly accessi- 
ble in experiments by using, respectively, time-of-flight 
techniques and two-photon Bragg spectroscopy 8]. The 
study of collective modes also provides a useful experi- 
mental technique to investigate the role of interactions 
and beyond mean-field effects 6] . Within a local density 
approximation (LDA) for systems in harmonic traps we 
calculate the frequency of the lowest compressional mode 
as a function of the interaction strength in the crossover 
from the TG gas to the super- Tonks regime. 

We consider a ID system of N spinless bosons de- 
scribed by the following contact-interaction Hamiltonian 
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where m is the mass of the particles, Zij — Zi~ Zj denotes 
the interparticle distance between particle i and j and 
giD is the coupling constant which we take large and 
negative. The study of the scattering problem of two 
particles in tight waveguides yields the following result 
for the effective ID coupling constant g\D in terms of 
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the 3D s-wave scattering length a^o 10] 
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where a± = ^Jh/muj^ is the characteristic length of the 
transverse harmonic confinement producing the waveg- 
uide and A = |C(l/2)|v / 2 = 1.0326, with ((■) the Rie- 
mann zeta-function. The confinement induced resonance 
is located at the critical value a 



a_\_/A and corre- 
sponds to the abrupt change of gm from large positive 
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The renormalization @ of the effective ID coupling con- 
stant has been recently confirmed in a many-body calcu- 
lation of Bose gases in highly elongated harmonic traps 
using quantum Monte Carlo techniques [Tl| . 

For positive g\D, the Hamiltonian JQl corresponds to 
the Lieb-Liniger (LL) model. The ground state and ex- 
cited states of the LL Hamiltonian have been studied 
in detail [12j and, in particular, the TG regime corre- 
sponds to the limit gm = +oo. The ground state of 
the Hamiltonian with gm < has been investigated 
by McGuire and one finds a soliton-like state with 
energy E/N = -mgl D (N 2 - l)/2Ah 2 . The lowest-lying 
gas-like state of the Hamiltonian with gm < corre- 
sponds to a highly-excited state that is stable if the gas 
parameter nam *C 1, where n is the density and a\o 
is the ID effective scattering length defined in Eq. J5J. 
This state can be realized in tight waveguides by crossing 
adiabatically the confinement induced resonance. The 
stability of the gas-like state can be understood from a 
simple estimate of the energy per particle. For a contact 
potential the interaction energy E int /N = gmng 2 (0)/2 
can be written in terms of the local two-body correlation 
function 32(0) = (ip* (z)^ (z)ijj(z)ip(z)) /n 2 , where ip>, ip 
are the creation and annihilation particle operators. In 
the limit gm — > —00 one can use for the correlation func- 
tion the result in the TG regime Q 32(0) ~ -k 2 n 2 a\ D I 'i , 
which does not depend on the sign of gm- In the same 
limit the kinetic energy can be estimated by Ek in /N ~ 
n 2 h 2 n 2 /(6m), corresponding to the energy per particle 
of a TG gas. For the total energy E = E^ iri + E int one 
finds the result E/N ~ Tt 2 h 2 n 2 / (<om) -TT 2 n 2 n 3 a 1D /(3m), 
holding for nam <C 1. For nam < 0.25 this equation of 
state yields a positive compressibility mc 2 = ndfi/dn, 
where u = dE/dN is the chemical potential and c is the 
speed of sound, corresponding to a stable gas-like phase. 
We will show that a more precise estimate gives that 
the gas-like state is stable against cluster formation for 
nam 0.35. 

The analysis of the gas-like equation of state is car- 
ried out using the VMC technique. The trial wave- 
function employed in the calculation is of the form 
iPt(z\, zjv) = Yii<j f( z ij)> where the two-body Jas- 



trow term is chosen as 
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The cut-off length Z is a variational parameter, while 
the wave vector k for a given Z is chosen such that the 
boundary condition at z = imposed by the 5-function 
potential is satisfied: — fctan(fcZ) = 1/am. For distances 
smaller than the cut-off length, \z\ < Z, the above wave 
function corresponds to the exact solution with positive 
energy of the two-body problem with the interaction po- 
tential gm5(z). For gm < (am > 0) the wavefunction 
f(z) changes sign at a nodal point which, for Z 3> am, is 
located at \z\ = am- The variational energy is calculated 
through the expectation value of the Hamiltonian (Q on 
the trial wave function E = (iI>t\H\iPt) / ('4 > t\' 1 I j t) ■ In the 
calculations we have used TV = 100 particles with peri- 
odic boundary conditions and because of the negligible 
dependence of the variational energy on the parameter Z 
we have used in all simulations the value Z = L/2, where 
L is the size of the simulation box. Calculations carried 
out with larger values of N have shown negligible finite 
size effects. 

The results for the variational energy as a function of 
the gas parameter nam are shown in Fig. ^ with solid 
symbols. For small values of the gas parameter our varia- 
tional results agree very well with the equation of state of 
a gas of hard-rods (HR) of size am (thick dashed line). 
The HR energy per particle can be calculated exactly 
from the energy of a TG gas by accounting for the ex- 
cluded volume @ 
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For larger values of na\o, the variational energy increases 
with the gas parameter more slowly than in the HR case 
and deviations are clearly visible. By fitting a polynomial 
function to our variational results we obtain the best fit 
shown in Fig. 2] as a thick solid line. The compressibil- 
ity obtained from the best fit is shown in Fig. ^ as a 
thin solid line and compared with mc 2 of a HR gas (thin 
dashed line) . As a function of the gas parameter the com- 
pressibility shows a maximum and then drops abruptly 
to zero. The vanishing of the compressibility implies that 
the system is mechanically unstable against cluster for- 
mation. Our variational estimate yields nam — 0.35 for 
the critical value of the density where the instability ap- 
pears. This value coincides with the critical density for 
collapse calculated in the center of the trap for harmoni- 
cally confined systems ^l|- It is worth noticing that the 
VMC estimate of the energy of the system can be ex- 
tended beyond the instability point, as shown in Fig. ^ 
This is possible since the finite size of the simulation box 
hinders the long-range density fluctuations that would 
break the homogeneity of the gas. This feature is anal- 
ogous to the one observed in the quantum Monte-Carlo 
characterization of the spinodal point in liquid 4 He [15J . 

As shown in Fig.^ the HR model describes accurately 
the equation of state for small values of the gas param- 
eter. A similar accuracy is therefore expected for the 
correlation functions of the system. The correlation func- 
tions of a HR gas of size am can be calculated from the 
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FIG. 1: Energy per particle and compressibility as a function 
of the gas parameter naio- Solid symbols and thick solid 
line: VMC results and polynomial best fit; thick dashed line: 
HR equation of state [Eq. Q]. Thin solid and dashed line: 
compressibility from the best fit to the variational equation 
of state and from the HR equation of state respectively. 



exact wave function ipHR — lli<j I ^[^(a' — Z j)/L]\, 
where the set of coordinates {z'j} is obtained from the set 
{zj} with the ordering Z\ < z 2 — oid < z$ — 2am < ... < 
zn — (N— l)a\D using the transformation z'j = Zj — jam), 
with j — 1,2, ...,7V. We calculate the static structure 
factor S(k), which in terms of the density fluctuation op- 



erator pk =J2i=i e%kZi > 1S defined as 
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and the one-body density matrix 
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Contrarily to the TG case, it is not possible to ob- 
tain analytical expressions for gi(z) and S(k) in the HR 
problem. We have calculated them using configurations 
generated by a Monte Carlo simulation according to the 
exact probability distribution function \iPhr\ 2 - The re- 
sults for the static structure factor are shown in Fig. [21 
Compared to S(k) in the TG regime, a clear peak is vis- 
ible for values of k of the order of twice the Fermi wave 
vector kp = nn and the peak is more pronounced as 
nam increases. The change of slope for small values of 
k reflects the increase of the speed of sound c with nam- 
The long-range behavior of gi(z) can be obtained from 
the hydrodynamic theory of low-energy excitations • 
For \z\ » £, where £ = Ti/{^/2mc) is the healing length 
of the system, one finds the following power-law decay 
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where the exponent a is given by a - 
TG gas mc = ir?in, and thus «tg — 
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= mc/(2irTin). For a 
1/2. For a HR gas 
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FIG. 2: Static structure factor S(k) for a gas of HR at differ- 
ent values of the gas parameter nam (symbols) and for a TG 
gas (dashed line). 
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FIG. 3: One-body density matrix gi(z) for a gas of HR at 
different values of the gas parameter nam (symbols) and for 
a TG gas (dashed line). 



one finds a = cvtg/(1 — nam) 2 and thus a > oltg- This 
behavior at long range is clearly shown in Fig. [3] where 
we compare gi(z) of a gas of HR with nam — 0.1, 0.2 
and 0.3 to the result of a TG gas 01 ■ The long-range 
power-law decay of gi(z) is reflected in the infrared di- 
vergence of the momentum distribution n(k) oc l/|fc| 
holding for <§C l/£. A gas of HR exhibits a weaker 
infrared divergence compared to a TG gas. The corre- 
lation functions of a HR gas at nam — 0.1, 0.2 should 
accurately describe the physical situation of a Bose gas 
with large and negative gm- For nam = 0.3 we expect 
already some deviations from the HR model, as it is ev- 
ident from the equation of state in Fig. ^ which should 
broaden the peak in S(k) and decrease the slope of the 
power-law decay in gi(z) at large distances. The analy- 
sis of correlation functions clearly shows that the super- 
Tonks regime corresponds to a Luttinger liquid where 
short range correlations are significantly stronger than in 
the TG gas. 
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FIG. 4: Square of the lowest breathing mode frequency, u> , 
as a function of the coupling strength Na^ D /a^ for the LL 
Hamiltonian (gw > 0) and in the super-Tonks regime (giD < 
0). The dashed line is obtained from the HR expansion (see 
text). 



Another possible experimental signature of the super- 
Tonks regime can be provided by the study of collec- 
tive modes. To this aim, we calculate the frequency of 
the lowest compressional mode of a system of N parti- 
cles in a harmonic potential V ex t = Si=i muJ z z i We 
make use of LDA which allows us to calculate the chem- 
ical potential of the inhomogeneous system /2 and the 
density profile n(z) from the local equilibrium equation 
fi = n[n{z)] + muj 2 z 2 Jl, and the normalization condition 

N = J R R n(z)dz, where R = ^2fl/(muj 2 ) is the size of 
cloud. For densities n smaller than the critical density 
for cluster formation, fi[n] is the equation of state of the 
homogeneous system derived from the fit to the VMC en- 
ergies (Fig. ^) . From the knowledge of the density profile 
n(z) one can obtain the mean square radius of the cloud 
(z 2 ) = J R n(z)z 2 dz/N and thus, making use of the result 



= -2 



d(z 2 )/dw 2 z 
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one can calculate the frequency u> of the lowest breathing 
mode. Within LDA, the result will depend only on the 
dimensionless parameter Na 2 D /a 2 , where a z — ^Jh/muj z 
is the harmonic oscillator length. For guj > 0, i.e. in 
the case of the LL Hamiltonian, the frequency of the low- 
est compressional mode increases from uo — \fiui z in the 
weak-coupling mean- field regime (Na\ D /a 2 3> 1) tow = 
2oj z in the strong-coupling TG regime {Na\ D /a 2 < 1). 
The results for to in the super-Tonks regime are shown 
in Fig. 21 as a function of the coupling strength. In the 
regime Na 2 D /a 2 -C 1, where the HR model is appro- 
priate, we can calculate analytically the first correction 
to the frequency of a TG gas. One finds the result 
u = 2uj z [1 + (16V^/15tt 2 ) (Na\ D /al) 1 ' 2 + •••]■ Fi S- H 
shows that this expansion accurately describes the fre- 
quency of the breathing mode when Na 2 D /a 2 <C 1, 
for larger values of the coupling strength the frequency 
reaches a maximum and drops to zero at Na 2 D /a 2 ~ 0.6. 
The observation of a breathing mode with a frequency 
larger than 2u z would be a clear signature of the super- 
Tonks regime. 

In conclusion we have pointed out the existence of a 
strongly correlated regime in quasi-lD Bose gases beyond 
the Tonks-Giradeau regime. This regime can be entered 
by exploiting a confinement induced resonance of the ef- 
fective ID scattering amplitude. We calculate the equa- 
tion of state of the gas in the super-Tonks regime using 
VMC and we estimate the critical density for the on- 
set of instability against cluster formation. The static 
structure factor and one-body density matrix are calcu- 
lated exactly within the hard-rod model, which provides 
the correct description of the system for small values of 
the gas parameter. For harmonically trapped systems 
we provide explicit predictions for the frequency of the 
lowest compressional mode. 
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